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Abstract 

It is shown that gauge theories are most naturally studied via a polar 
decomposition of the field variable. Gauge transformations may be viewed 
as those that leave the density invariant but change the phase variable by 
additive amounts. The path integral approach is used to compute the par- 
tition function. When gauge fields are included, the constraint brought 
about by gauge invariance simply means an appropriate linear combina- 
tion of the gradients of the phase variable and the gauge field is invariant. 
No gauge fixing is needed in this approach that is closest to the spirit of 
the gauge principle. We derive an exact formula for the condensate frac- 
tion and in case it is zero, an exact formula for the anomalous exponent. 
We also derive a formula for the vortex strength which involves computing 
radiation corrections. 

1 Introduction 

The density phase transformation for bosons is quite well-known to those who 
work in the fields of superfluidity [1] [2] and other related areas such as supercon- 
ductivity and possibly even quantum optics. Jackiw and collaborators [3] have 
recently introduced some of these ideas independently in the context of rela- 
tivisitic quarks. Since the subject of bosons is vast, we refer the reader to the 
review by Ceperley[4] for further references. However a few basic works do de- 
serve mention. First there is the work by Penrose and Onsager[5] on estimating 
the superfluid fraction in He'^. Then there is the Hohenberg-Mermin- Wagner 
theorem[6] that precludes long-range order with a broken continuous symmetry 
in two dimensions. Finally, there is work by Ceperley[4] using Quantum Monte 
Carlo. 

The main point of this article is to highlight the usefuUness of the hydro- 
dynamical approach in tackling gauge theories. We argue that this approach is 
the natural way in which to study gauge theories. Although by itself none of 
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the results in this article are particularly new, the novel aspect is the manner 
in which the theory is formulated that renders it susceptible to generalization 
in future publications to fermions coupled to gauge fields and later on to rel- 
ativistic gauge theories and finally to nonabelian relativistic gauge theories. 
Much of this agenda has already been completed by the high energy community 
recently, specifically, the work of Jackiw and Polychronakos [3] in introducing 
anticommuting Grassmann variables into a fluid dynamical description will be 
very relevent to our future work. But our thrust is to compute microscopic 
correlation functions rather than just develop formalism. 

One concrete result in this article is an exact formula for the condensate frac- 
tion. Since the condensate fraction depends only on the asymptotic properties 
of the one-particle Green function, and this is given exactly via bosonization, we 
may conclude that the derived result for this quantity is therefore, exact. The 
formula for the condensate fraction may also be derived using other methods, 
most recently Liu [7] has used the eigenfunctional theory to derive an expression 
for the same quantity that agrees with the expression derived in this article. 

We apply this method to compute the velocity-velocity correlation function 
and demonstrate the power of this method in accounting for vortices in a gauge 
invariant manner. 



2 Density Phase Transformation 

Consider the action of nonrelativistic spinless bosons. We use units such that 
h = 2m=l. 

Sfree = ^ ^ dt J d'^X V^^x) {iOt + V^) ^(x) (1) 

In the path integral approach, ?/'(x) is just a complex number that is defined at 
each point x. Every complex number can be decomposed into a magnitude and 
a phase. Using the ideas explained in our earlier work [8] 

V(x) = e-' (2) 
Prom this we may write down the current, 

J(x) = -p(x) vn(x) (3) 

-i(3 



Sfree = J dt J d'^'X 



(4) 



pdtii-p{viif- — 

In deriving this action we have already used some boundary conditions. From 
the text by Kadanoff and Baym[9] we learn that the Green functions of bosons 
obey the KMS (Kubo-Martin-Schwinger) boundary contitions. These boundary 
conditions translate in the path integral representation to, 

V'(x,i-i/3) = 6'^" V(x,t) (5) 
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This in turn means that the number conserving product is invariant under this 
discrete time translation. To see this we examine, 

Observe that V''' (x, i - = [ip{-x,t + i/S)]'^ = e''^'' ^'^(xjf) and -ipix ,t - = 
e^'^ i^{x,t). Thus, V^(x,f - il3)il){yi ,t - if3) = i/;^ {x,t)^{x ,t). As pointed 
out in an earher work[8] the phase variable may be written as a sum of two 
terms a position independent term which is the conjugate to the total number 
and a position dependent term that is related to currents and densitites. Thus, 
n(x) = + n(x). We have just argued that n(x) and p(x) are invariant under 
the discrete time translation. Thus in order to preserve the KMS boimdary 
condition we must impose Xo{t — ij3) = Xo(i) + i /3/x. Therefore we find that 
the conjugate to the total number makes its presence felt in a very nontrivial 
manner. The boundary condition that has been used in deriving the action is 
N{—ifi) = iV(0) where N is the total number of particles. 



3 Gauge Transformations 

Here we examine what sorts of changes are brought about by the imposition of 
local gauge invariance. Gauge transformations in the usual language is given 
by, 

V. (x) = ^ ^ V(x) (6) 

\ = Ao + dtA (7) 

a' = A + VA (8) 

In order to find an action invariant under these transformations, we replace 
derivatives by covariant derivatives (minimal coupling). Thus we have — > 
D/j, = — ie A^. In the usual language, 

S = J ^^idt + eAo + {V -ie Af) V (9) 

Here and henceforth by / we mean Jg dt J d^x Now we rewrite this as a 

sum of two parts the free term plus the term that couples to the gauge fields 
and finally the term involving only gauge fields namely, the curvature term. 

S = Sfree "I" ^int + Sgauge (l^) 



^free — 



pdtii- p(vn)2 - 



2_ (M 

4p 



(11) 



S^nt^ j peAo-e^ j pA^-2e J p (A ■ VH) (12) 
= -lJ (13) 



Sgauge 



This action is gauge invariant provided we have 

n'=n-eA; p = p (14) 
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4 The Jacobian Determinant 



In this section we attempt to rewrite the partition function in terms of the 
density and phase variables. Consider the partition function in the original 
Bose language. 

Z= I D[il)] D[i,^] e'^!'''''''^! (15) 



(16) 



In the density phase variable language we have, 

Z = j D\p\ D[n] J(n,p) e^^l"'") 



Here J(n, p) is the appropriate Jacobian determinant that tells us how the 
measure transforms. Fortunately, this Jacobian is a constant. To see this we 
write the definition of J as, 



j(n,p) 



sn 



Sp 



(17) 



In our next article, wc shall sec that the density phase variable ansatz (DPVA) 
for fermions[8] which includes a phase functional that is nonlocal in position 
space also leads to a Jacobian that is constant. This can be checked using the 
Mathematica software. 



5 Propagator With Two-Body Forces 

Here wc compute the free propagator using path integrals in the density phase 
variable language. This is done to convince ourselves of the basic soundness of 
the approach. Also we operate in the limit where the mean density is constant 
and so that we may ignore the density fluctuations in the long- wavelength limit. 
This will be made precise soon. The main advantage of this approach is the 
ease with which we may treat interactions. Just to highlight this fact we also 
include density-density interactions. 



Sfr 



(18) 



4po 

Using the boundary conditions we may write, 



V 



H(x, t) = -nt + Y, e'^-^X^n e^"* (20) 
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The action then may be written as foUows. We have added a two-body potential 
since we may do so without additional effort. 

SfuU = i Pl^ PO,0 + ^(-^P Zn)PqnXqn 
q,n 

+i/3iV° J2 ^q,n^-q,-n + Pq,nP-q-n + |^Pq,nP-q -n 

q,n q,n q,n 

(21) 

Here Zn — 21171/ p. From this we may write down a formal expression for the 
propagator. We set p ~ po in the argument of the square roots. 

G{x,t) = -ipo e'^^exp 



Here Wq = q^(q^ + 2poVq) is the Bogoliubov dispersion. The Matsubara sums 
may be performed quite easily and we may derive an expression for the equal- 
time version of the time-ordered propagator, 

^0 + S<i., 0) = E (^(S!±|2!^ _ 1 j (1 _ eo.(q.x)) (23) 

where, 

G< (x, 0) = e-'"+'''^^ G< (x, 0) (24) 

The main reason why the one-particle propagator is interesting is because while 
the Gaussian approximation leads to just Bogoliubov's theory as far as the 
the computation of density-density correlation functions arc concerned, the one 
particle properties are singular in one dimension, just as in the case of fermions 
in one dimension. 



(2 - e" 



q,n 



(22) 



5.1 The Condensate Fraction at Zero Momentum 

The analysis in the preceeding sections employs a Gaussian approximation which 
is valid if the density fluctuations are small compared to the mean density. To 
see this more clearly we write the condition as follows, 

V< PqP-q > « (25) 

where is the total number of particles. The density- density correlation is 
given by < PqP-q >= S{q}. Using the same Gaussian approximation we 
may deduce that 

(q^ + 2poVq) 2 
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Thus this scheme is self-consistent (as opposed to self-contradictory) only if, 



r < iV° (27) 

(q2 + 2/90 Wq) 2 

If we assume that = Alql"*"*"^ we have, the condition( for iV° ^ 1 ), 

(2poA) |qr » -1 (28) 

It would appear that so long as A > and m < — 2, this holds for small |q|. 

Thus the delta- fimction potential ( m = —2 ), the Coulomb potential in 2d ( 
m = — 3 ) and 3d ( m = — 4 ) all obey the inequality for small |q|. For large 
enough |q| the left hand side is zero but >• — 1 is not acceptable, hence the 
Gaussian approximation breaks down for large |q|. Thus it would appear that 
the results for the correlation functions are exact in the asymptotic limit. In real 
space, this is the |x| oo limit. Let us examine this limit of the one-particle 
Green function. To this end let us write a formal expression for the momenum 
distribution as 

nk = /o N° 5k,o + /k (29) 

here /k is a continuous function of |k| and is of order unity. < /o < 1 is the 
condensate fraction and N'^ is the total number of particles. The propagator in 
real space is then given by, 

< V^(x, 0)V(0, 0) >=foPo + ^Yl cos(k.x)/k (30) 

k 

In the ultra-asymptotic limit |x| — > oo, the cosine may be set equal to zero and 
we have just the first term. Thus it would appear that the Gaussian approxi- 
mation being exact in this limit, yields the exact condensate fraction but not 
the exact /k. This is analogous to the assertion that in the Fermi case, the 
analogous method yields the exact quasiparticle-residue but not the full mo- 
mentum distribution and in case the quasiparticle-residue is zero it yields the 
exact anomalous exponent. Here too when the condensate fraction is zero, we 
have to instead compute the anomalous exponent. Thus the exact condensate 
fraction is given by. 



/o = exp 



4iV0 I |q| 



(31) 



For the delta-function interaction in Id, the condensate fraction is zero. This 
could not have been guessed from Bogoliubov's theory that predicts (rather 
assumes) that a condensate always exists. In two dimensions, for the gauge 
potential Vg = const. /q^, the condensate fraction is zero. We could now evaluate 
the total energy of the system (per particle) and compare with the results of 
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Lieb and Liniger [10] who solve their model in Id or with the results of Schick[ll] 
which is valid in 2d but we do not expect the comparisons to be favorable since 
for the total energy to come out right wc need the exact /k rather than /o which 
does not contribute at all. Since the Gaussian approximation does not give us 
this we shall not bother performing this calculation. 

5.2 Anomalous Exponents 

In case the condensate fraction is zero, we have the bosonic analog of the Lut- 
tinger liquid or the Lieb-Liniger liquid. Here we have to instead compute the 
anomalous exponents which are also given exactly in the Gaussian approxi- 
mation. Interestingly we may address the question whether the condensate is 
destroyed in more than one dimension also. In passing we note that the main 
purpose of this article is to hint at the usefullncss of this approach in studying 
fermions coupled to gauge fields. There we have to polar decompose Grass- 
mann variables - an exercise which is still in its infancy, although considerbale 
inroads have been made by the author. For a delta- function interaction in Id it 
is well-known that there is no condensation and instead we have to compute the 
anomalous exponent of the propagator. It is clear from the preceeding sections 
that we may write, 

7o + 5<(x,0)«--^ r dq^^^^f^il-cos{q.x))^-jLn{\x\) (32) 
47rpo Jo \q\ 

where 7 = {2vopQ)i /{Airpo) is the anomalous exponent. The interesting ques- 
tion is whether wc can have the destruction of the condensate in two or three 
dimensions for realistic potentials ? This would be the bosonic analog of the 
question 'Does Fermi liquid theory break down in two or three dimensions ?' 
In two dimensions, for the interaction Vq = const. /q^ (gauge potential) the 
condensate fraction is zero. It would appear that for realistic Fourier trans- 
formable potentials in three dimensions we have Bose condensation but not in 
less than three dimensions. We have derived an exact formula for the conden- 
sate fraction that is valid for functional forms of the interaction that are Fourier 
transformable and are subject to the constraint mentioned previously namely 
that the long-wavelength limit is exactly given by the Gaussian approximation. 

6 Correlation Functions With Gauge Fields 

When one is performing the path integral with gauge fields, one must be care- 
ful about preserving gauge symmetry. This is the crucial aspect that may be 
elegantly treated in the present approach. The usual method for treating such 
constraints in the context of path integrals is the well-known Faddeev-Popov 
method [12]. However we are able to treat gauge symmetry in all its generality 
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without ever having to fix the gauge at any time. Thus the path integral is to 
be performed such that the following 

e ^0 + 9t^' =eAo + dtU (33) 

e a' + Vn' = e A + Vn (34) 

constraints are obeyed. We may therefore simply solve for the gauge fields 
straight away in terms of the conjugate variables. Define an arbitrary gauge 
constant C. This is constant in the sense that changes in the phase of the field 
cancel out the changes in the vector potential. However, in the end we will have 
to integrate over this variable as well. It will be shown that in the case when 
e = we recover the noninteracting propagator. 

Co = eAo + d,U (35) 

C = e A + vn (36) 

In the noninteracting limit (e = 0), 11 and 11 differ at most by a trivial constant 
and hence the vector C is irrotational. This observation will be made use of later. 
The action now may be recast in terms of the gauge constant and the density and 
phase variables. The field tensor is then simply given by, F^i, = d^Ci, — d,^C^. 
Thus the total action simply reads, 

We may see by examining Eq.( 37) that in the limit e ^ 0, the field tensor has 
to be zero for the partition function to be nonvanishing, thus in this limit we 
must have =- d^jj\ for some scalar 11 so that F^^v = 0. This then gives us the 
action functional of the free theory as it should. 

The next question that is worth answering is the following. How should 
the current operator be defined when there are gauge fields ? In particular is 
the current operator gauge invariant ? The density surely is gauge invariant. 
If we would like to treat currents and densities as part of a canonical set of 
variables then we better have a definition that is gauge invariant. The definition 
J = — pVn is not gauge invariant. Indeed, it changes by additive amounts 
everytime a gauge transformation is performed. In order to remedy this we 
redefine the current to be J = —p C. As we have seen, in the noninteracting 
limit, this reduces to the form already shown, since there C = VH. In general 
however, it does not. Now we would like an expression for the field operator 
also in terms of C. Using the DPVA[8] in terms of currents and densities we 
may write, 

V'(x) = e-'^°+'r*^VpW (38) 
Using the formula for the current in terms of the gauge constant we may write, 

V;(x,i) = e-'^°(*^-'r*''("''*'v/p(^ (39) 
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Computing the propagator of the interacting theory then means evaluating the 
following path integral. 

G(x, i) = 9 / ^[^]^ W ' ^(^' ^)^^ (0' 0) (40) 

and 

Z = j D[C]D[p] e' ' (41) 

where S is given as in Eq.( 37). The integration over the gauge constant C may 
be justified as follows. In the original path integral we had to integrate over tp, 
and the four 's. Thus we had six variables in all, except that the gauge 
constraint meant that one of the variables, say the phase of the field depended 
on the gauge fields thus the number of independent variables reduced to five. 
In the definition Eq.( 41) we have a similar situation - integrating over the four 
components of the gauge constant plus the density p makes five variables as it 
should. The action of the free gauge field is given by, 

S„ = ^J (doC - VCf + ±JcV'C+^JiV-Cf (42) 

A word regarding units is appropriate. We have set h = 2m = 1. The presence 
of gauge fields implies that there is a further dimensionful parameter namely 
the speed of light which wc denote by vi . If we set vi = I then all quantities are 
dimensionless. This is undesirable. Henece we retain the speed of light as it is. 
The various quantites in this work now have the following dimensions, [pkn] = 1 
but [p(x,t)] = ; [C(x,t)] = [Ck„] = [k] = [V] - [vi] ; k^ = [Energy] = 

[do] ; (3 = [Energy]~^ and finally [e'^/V] = [Energy]'^. Using the decomposition 
into the various modes we may write, 

^ = E Pk,„P-k,-„ + i-ip) ^ Pk,„CO,„ + (iPN^) Ck,„ ■ C_k,-„ 

k k,n k,n 

+(^^)^ E 4Ck,n-C_k,-„-(*/3)^ ^(2zz„ k)Ck,„C°k,-„-(i/3)^ E k'Ck,nC-i 

k,n k,n k 

+ ^E'^' Ck,„-C_k,-„-(z/3) vf ^^(k-Ck,„)(k.C_k,-„) (43) 

k,n k,n 

After integrating out the (7° we have, 

^ = E + ^k^) ^k,nP-k,-„ + (iP) ^Pk„(k • Ckn) 

kn ^ ^ k,n 

+{iP)Y(N°+vf^kA Ck,„-C_k,-„-(i/3)i;f ^l](k-Ck„) (k-C_k,-„) 

kn kn 
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+ E • C_k,-„ - E § ■ Ckn)(k • C_k,-„) (44) 

kn kn 

Observe that in the hmit e — > 0, we must have Ck,n = k Xk,n for the partition 
function to be nonvanishing. Thus wc recover the nonintcracting theory in this 
limit. Now we would like to compute two quantities. One is the dynamical 
density-density correlation function. Here we would like to see radiation cor- 
rections to corresponding spectral function. However, it is not present at the 
Gaussian level at which we are presently operating. Thus we shall be content 
at reproducing results equivalent to Bogoliubov theory. This is a novel way 
of recovering the Bogoliubov spectrum by introducing gauge fields and using 
the path integral formalism. The other quantity of interest is the circulation of 
the gauge constant. This quantity is zero in the nonintearcting case since the 
vector C may be expressed as the gradient of a scalar. In fact we may see from 
the expression of the current that the velocity of the fluid is simply given by 
V = — C. This is irrotational in the absence of gauge fields. In general however 
this is not the case. The circulation of the velocity around a closed loop P is 
given by, V{P, t) = — fpdl- C(x, t). Physically, the quantity V{P, t) is a mea- 
sure of the strength of vortices in the system. The average < V{P, t) >= is 
trivially zero. Thus wc have to examine the fluctuation < V'^{P,t) >2. First 
the density-density correlation function. In order to evaluate this we must first 
perform the integral : J D[C] e'"^. Again here we may use the usual procedure 
of translating the integration variable by a constant amount and we obtain the 
following effective action. 

S=Ji,Vfk^ + imS^V,»P-.,-» (45) 



47VO ^ V k2 

k.n ' 

From this we may immediately deduce the static structure factor as 

5(k) = - (46) 
where the Bogoliubov dispersion Wk is given by the positive real solution to, 

k2 + ^^^^^±^=0 (47) 

In other words, 

(^k= |k|v/k2 + 2po«k (48) 

where Wk = e^/k^. Therefore as is well-known, independent of the dimension- 
ality of space, the interaction in the presence of gauge fields is forced to be 
of the form ~ l/k^. In three space dimensions, this corresponds to Coulomb 
interaction. This is a necessary consistency check. 

Some may want to evaluate the propagator. For this we have to evaluate the 
correlation function < C'(x, t)C^{0, 0) >. There is a more interesting reason to 
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study this, namely to compute the vortex strength. We have found that even 
though at the Gaussian level the density-density correlation is unremarkable, 
the velocity-velocity correlation function does exhibit some new physics. To see 
this we note that integrating out the velocity variable means replacing Ckn — *■ 
Ckn + kAkra for an appropriate A. This makes all the terms involving the speed 
of light and most other terms drop out and we are led to Eq.( 45). However if 
we integrate out the density first and retain the velocity as it is then we find 
that the final action is no longer so simple. In particular, it will involve both 
longitudinal and transverse terms, the latter being responsible for vortices as 
we shall see. If we integrate out the p first we are led to the following effective 
action. Here we have also introduced a source for velocity. 

S, = -z/3iV° ^ (k • Ck„)(k . C_k,-„) + Ak„ • Ck„ 

+{i(3)Y^(N^ + v^^\A Ck,„-C_k,-„-(i/3)^;f ^Y.^^'^^^^) (k-C_k,-„) 

kn kn 

E 'I Ck„ • C_k,-„ - E § ■ Ckn)(k • C_k,-„) (49) 

kn kn 

It can bee seen from Eq.( 49) that in the limit t;; ^ oo, the partition function 
is nonvanishing only if k^Ckn • C_k,-n = (k • Ckn)(k • C_k,-n)- In other 
words, only if Ckn = kAkn- This is the same as saying that the velocity is 
the gradient of some scalar v = — VA. Thus in this limit there are no vortices. 
Vortices arise as a result of a finite speed of light, in other words due to radiation 
corrections. Now we evaluate the partition function 2'([A]) = Z([A])/Z([0]) 
where ^([A]) = / D[C\ e*'^°. After some tedious algebra we arrive at. 



Z{A) = exp 
where, 



\Yf{Kn){k ■ Ak,„)(k- A_k,-n) - i^5(k,n)Ak,n • A_k-„ 

(50) ■ 



kn 



kn 



/(k,n) = (2/3ArO)-i 

7-0\-l 



(iz„/k2)2 1 _ z2 1 



k2 + 2p0e2/k2 ^ 2p0e2^ ' k2^ 



4/k^ 



k2 + 2p0e2/k2 



g{k,n) = (2/3Ar")-^ ( 1 + ^(^fk^ + <) 



(51) 
(52) 



11 



We may deduce a formula for the time ordered vefocity-vefocity correlation 
function, 



<^T C*(x ,0 CJ(x,t)^ = [V\Vi F(x -x,t -t)+H{x -x,t - <5i,jj 

, , (53) 

F{x',t') = ^/(k,n) e'l'-^'e-^"*' ; H{x' ,t') = ^(?(k,n) e^'^-^'e-^"*' (54) 

k.n k,n 

The S and below are defined in the appendix. 

S(t -t;P, P) = -i i f di ■ dl H{x' -x,t - t) (55) 

Also we may write for the propensity to create vortices (defined in the appendix), 

^2 r r . p-A|x'-x| 



AlTV 



(t(P,P) = -3^ f i dl-dl— (56) 



where A = {2p^e^ /vf)^ . This quantity looks finite. Hence there is no vortex 
instability in the Gaussian approximation. Let now compute the vortex strength 
as defined earlier. 

{V^iP, t))= i i di ■ dl H{x' - X, 0) (57) 
J\'eP Ji€P 

k 

where a';(k) = (2p°e^ + vfls.'^)i . This dispersion shows that the photons have 
acquired mass by coupling to matter. An exphcit evaluation of this and further 
analysis is possible and will not be done here since our interest is merely to 
highlight the usefulness of these ideas. However some intutive justification of 
the above formulas is in order. First the vortex strength and suceptibility are 
zero for a noninteracting system. This is hardly surprising since the velocity 
operator of free bosons is irrotational in our formalism. The author has tried 
to prove this independently using the second quantization formalism but has 
been unsuccesful in proving it rigorously. It appears that the only obstruction 
to this conclusion at least at the formal mathematical level, is the lack of self- 
adjointness of the canonical conjugate to the total number of bosons [8] [13]. In 
the high density limit this is not a problem. Thus as far as the asymptotics are 
concerned, our approach provides reliable answers. In other words, the expres- 
sion for the velocity-velocity correlation function is valid in the largo scperation 
|x — X I ^ cxD limit. However the main message is that the vortices are due 
to the finiteness of the speed of light. That is, even in an interacting system, 
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there are no vortices if we include only the longitudinal density-density inter- 
actions. Coupling to transverse radiation fields is needed to generate vortices. 
This may be understood by realizing that in order to create vortices we have to 
supply angular momentum. Photons being spin-1 particles, possess the angular 
momentum that the charged bosons use to generate vortices. 

7 Conclusions 

We have written down formulas for the one-particle Green function of a ho- 
mogeneous Bosc system that is exact in the asymptotic limit. From the exact 
asymptotic form we may extract an exact formula for the condensate fraction 
and if it is zero an exact formula for the anomalous exponent. We have also com- 
puted the vortex strength and shown that radiation corrections arc responsible 
for the vortices. No other quantity is given exactly in the formalism outlined. 
The total energy per particle, roton minimum and other important physical 
attributes will have to involve going beyond the Gaussian approximation. Nev- 
ertheless, the agenda in the immediate future is to apply these ideas to fermions 
and then write down a theory of neutral matter with nuclei and electrons treated 
on an equal footing. 
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In this section, we derive a Kubo-like formula that relates the propensity for 
vortices emerging ( an analog of d.c. conductivity ) with microscopic velocity- 
velocity correlation functions. Let us write down the following formula for the 
vortex strength with a source for vortices in the interaction representation. As 
defined in the main text, the vorex strength may be written as follows. 



9 Appendix 




(59) 





(60) 



(T S) 




(61) 
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Here W{P, t) is a whirlpool source for the vortex V{P, t). From linear response 
theory, we may expect that for weak sources the vortex strength is proportional 
to the convolution of the source with a linear response coefficient. 

{V{P,t))=Y, dtY.{t-t-P,P)W{P ,t) (62) 
p, Jo 

Prom this we see that, 



Y.{t-t:P,P) = -i (T V{P , t )V{P, t)j (64) 
If the whirlpool source is time-independent, we may write, 

{V{P))=Y,'^iP,P')W{P') (65) 
p' 

a{P,P')= dtY:{t-t;P,P') (66) 



If cr(P, P ) = oo then it signals instability to vortices, namely arbitrarily weak 
pertubations can lead to vortices. Thus we have, 

T.{t-t;P,P') = -i j ^d\ dlj (t C\x' ,t)C^{x,t)'^ (67) 
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